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,i^ ' Abstract. We study the ellipticity probability distribution function (PDF) of voids in 

M-i. redshift space with galaxies as tracers of the shapes of voids. We find that the redshift space 

O ! distortion on the shape of voids statistically increases the ellipticities of voids, and leaves 

a prominent feature on the ellipticity PDF as a substantial reduction in the probability of 

j^ ■ having voids with small ellipticity. The location of this characteristic cutoff of the ellipticity 

PDF is an explicit function of the logarithmic growth rate, f{z) = ^"^J-^) ; and it can be 

used as a probe of cosmology once the radial density profile of voids is better understood. 

^ ' However, the biggest limiting factor for the use of ellipticity PDF as a probe of cosmology lies 

^is . in the Poisson noise from a small number of galaxies to define the shape of a given void. This 

00 ! Poisson noise creates a significant contamination of the resulting ellipticity PDF so that the 

^^ ■ shape of the original PDF is almost washed-out. Nevertheless, there is a way to overcome the 

cn . Poisson noise via the Alcock Paczynski test on the shape of stacked voids. In redshift space, 

since the void is elongated toward the line of sight, the stacked void has non-zero ellipticity, 

which can be a tell-tale of the logarithmic growth rate. Although some useful information 

of void ellipticity will be lost by stacking, in this way, we can see the effect of redshift space 

distortion as a source of anisotropy in the stacked void ellipticity. We think that the stacking 

analysis of the voids in redshift space is potentially a powerful tool to probe the cosmology. 
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1 Introduction 

The abundance of on-going and future galaxy surveys is pointing to the deeper understanding 
of the nature of accelerated expansion of the universe, as discovered via observations of 
luminosity distances to Type la supernovae [30, 32]. The main observables of galaxy surveys 
are characteristic length scales encoded in the matter power spectrum, P{k), such as the 
comoving Hubble horizon, H{z), and the angular diameter distance, Da{z) [37]. 

Not only can we measure geometric distances, DAiz) and H{z), from a galaxy survey, 
but we also measure the growth rate of the matter density fluctuation, f{z) = ^"^ t^f , via 
redshift space distortions [10]. 

The matter density fluctuation grows via gravitational instability competing against the 
expansion of the universe. One can obtain the growth rate by solving the following differential 
equation [18, 24, 40], 
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As we see in the eq. (1.1), growth rate has an explicit dependence on the dark energy equation 
of state. Therefore, by measuring the growth rate, one can obtain the information about the 
dark energy, thereby, the nature of accelerated expansion of the universe. 

Predictions for the growth rate from different theories of gravity and the nature of dark 
energy can be well parametrized in the simple form of f{z) = Vlm{zy ■ For example, 7 = 4/7 
for a ACDM model [11, 29], where the cosmological constant plays the role of observed 
accelerated expansion of the universe, and 7 = 0.68 for the DGP model of modifications 
of gravity [23]. Thus, it is crucial to measure f[z) from a given galaxy power spectrum to 
observationally unveil the nature of this mysterious component of the universe. Although a 
galaxy power spectrum gives a way to measure f{z) at different redshifts, due to the vast 
parameter space allowed, an accurate measurement of f{z) from a galaxy power spectrum 
alone is difficult [38]. Especially, when we use the redshift space galaxy power spectrum to 
the linear order, it is inevitable to have the growth rate, /(-z), be degenerate with a linear 
galaxy bias, 6l, such that Pg{k,fj,) = Pm{k,fi){l +/3^^)^, where /3(z) = f{z)/bL [16]. Thus, 
it is important to have an alternative method to measure the growth rate in addition to a 
redshift space galaxy power spectrum. Here, we propose the shape of voids in the redshift 
space as a sensitive probe of linear growth rate, f{z). 

The shape of the voids in redshift space has been studied in literature [25, 33, 34]. 
According to their study, the void volume increases in the redshift space as its boundary is 
stretched along the line of sight by redshift space distortion. This elongation of the void vol- 
ume in the redshift space is diametrically opposite to what would happen to the distribution 
of galaxies in overdense regions, due to a reversed sign of the equation of motion sourced by 
local density field [14] (i.e., 6 > for over dense regions and 6 <0 for voids). Schmidt et al. 
[35] studied the effect of redshift space distortion on the void probability function (VPF). 
The VPF is a measure of probability, Pq{V), that a randomly placed sphere of volume V 
contains no galaxy within [41]. For a given void in the real space, redshift space distortion 
stretches the void along the line of sight, increasing the volume of voids; thus, VPF, PQ(y), 
increases toward a larger V. 

Here, we study the ellipticity of voids in the redshift space as a probe of the linear 
growth rate, f{z). Specifically, we shall use the ellipticity probability distribution function 
of Park & Lee [27] (hereafter, PL07) and Lee &: Park [21]. PL07 proposed a void ellipticity 
distribution function as a sensitive probe of cosmology, noting that the shapes of voids are 
modulated by the competition between the tidal distortion and the cosmic expansion. Their 
model is based on the assumption that the underlying tidal field can be well described by the 
Zel'dovich approximation. They tested the model against the results from the Millennium 
simulations [39], finding a good agreement. 



In practice, defining a void from a distribution of galaxies is non-trivial, as the definition 
of a void differs from one void finding algorithm to another (Colberg et al. [4] and references 
there in). Regardless of the ambiguity in the exact definition and the boundary shape of a 
void, the method of PL07 gives a robust measurement of ellipticities of voids (i.e., defined 
as space around local density minima): it uses the distribution of galaxies inside a void as 
a tracer of the underlying shape of the tidal field. Void statistics have been studied both in 
simulations and galaxy redshift surveys [8, 12, 26, 28]. With the abundant data of the galaxy 
redshift survey available from on-going and future galaxy surveys, properties of voids can be 
better understood. Furthermore, voids found from these surveys will provide a unique way 
to probe the history of the structure formation, and henceforth, the nature of dark energy. 

Throughout this paper we use the maximum likelihood cosmological parameters of Ko- 
matsu et al. [18] (WMAP-hBAO-^SN of Table 1). 

The goal of this paper is to provide an analytic insight into the observed ellipticities of 
voids, and establish basis for their application to cosmology. In § 2, we briefly review notations 
and definitions used in PL07 for a real space void ellipticity probability distribution function. 
In § 3, we study the effect of linear galaxy bias on the void ellipticity PDF. In § 4, we study 
the effect of redshift space distortion on the shape and ellipticity of observed voids. In § 5, we 
study the effect of the Poisson noise on the void ellipticity PDF, which arises as a consequence 
of tracing a void shape by a limited number of field galaxies inside the void. We also show 
that this Poisson noise can be the biggest limitation on the use of PLOT method for voids 
with a small number of field galaxies inside. In § 6, we compare the analytic prescriptions 
derived in the previous sections against the void ellipticity PDF extracted from the galaxy 
catalog of the Millennium simulation. Finally, our conclusions are in § 7. 

2 Void in the real space 

In this section, we briefly review the method of PL07. First, voids are extracted from a given 
galaxy distribution using a given void finding algorithm. We then find the ellipticities of 
voids, and construct a void ellipticity distribution function. 

For each detected void, we define void galaxies, which reside inside a void, and calculate 
the inertia tensor of equal weight using A^vg void galaxies, 

Afvg 
a=l 

where Xi is the distance from the center of a void (i.e., the mean, Xj, has been removed). 
Suppose that lengths of the semiaxis of the best-fit ellipsoid to a void are pi, p2 and ps, 
ordered as Pi > P2 ^ Ps ^ 0. Then, one finds that they are proportional to the square root 
of the eigenvalues of the inertia tensor, /i, I2 and I3, ordered as Ii > I2 ^ I3 ^ 0. This 
motivates our defining the ellipticity of a void such that. 




V^l-- = l-\^. (2.2) 



For a fixed volume of an ellipsoid, V = -j-pij'2P3, we can exhaust all the possible shapes 
of the ellipsoid in terms of e and rj. Figure 1 shows five special configurations of ellipsoids 
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Figure 1. We show five special configurations of ellipsoids: (a) sphere, (b) filament, (c) pancake, (d) 
prolate spheroid and (e) oblate spheroid. In a plane of e — 77, the lower triangle shows the possible 
configurations of ellipsoid, where each corner corresponds to a special configuration such as a sphere 
(lower left), a filament (upper right) and a pancake (lower right). Genuine triaxial ellipsoids are 
located in the middle of the triangle. 



on the e-T] plane. Note that according to our definition (i.e., Pi > P2 ^ Ps ^ 0), we have 
< ?? < e < 1. 

(a) For e = r/ = (i.e., pi = P2 = Vi)-, we have a spherical void. 

(b) For e = r] = 1 (i.e., Pi ^ P2 = Ps = 0), we have a filamentary void. 

(c) For e = 1 and 77 = (i.e., Pi = P2 ^ Ps = 0), we have a pancake void. 

(d) For e = r] (i.e., pi > P2 = Ps), we have a prolate spheroid. 

(e) For r] = (i.e., pi = P2 > Ps), we have an oblate spheroid. 
We define the triaxiality parameter following Franx et al. [9], 

7?(2 - r?) 
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and further divide the triaxial ellipsoid into two categories: prolate ellipsoids (T > 0.5) and 
oblate ellipsoids (T < 0.5). 

2.1 Real-Space Ellipticity PDF 

Here, we calculate the ellipticity probability distribution function of PLOT based on the un- 
conditional joint probability density distribution of eigenvalues of a tidal tensor, {Ai, A2, A3} 

[7]. 



unconditional PDF is given by 
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where -fCi = Ai + A2 + A3 , ^2 = A1A2 + A2A3 + A1A3, and or^ is the rms density fluctuation 
smoothed at the Lagrangean void size, i?^, with a top-hat window function, W{kKi)^ 
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Here, we calculate A^(fc) = ^^2 with a linear Boltzmann code, CAMB [22], with the 
best-fit cosmological parameters of WMAP-FBAO-hSN [18]. 

We derive a Lagrangean void size, i?i, from conservation of the number density within 
the volume elements of Eulerian and Lagrangean space, 

n{x,z)(fx = n(fq, (2.6) 

where n{x,z) = n(l -|- 5{x,z)) is a real-space comoving number density of matter, n is a 
mean comoving number density, and x and q are Eulerian and Lagrangean coordinates, 
respectively. Therefore, we have 

i?L(q, z) = Re{x, z){l + 5(x, z))V3, (2.7) 

where Re is an Eulerian size of void. PL07 derived the conditional PDF of {Ai, A2} for a 
given density contrast, 6 = X^j=i Aj as follows, 
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Here, as studied by Lavaux & Wandelt [19], signs of eigenvalues, Aj tell us whether the void 
is spatially expanding or contracting along the corresponding semiaxes. In our notation, a 
given void spatially expands when Aj < 0, and contracts when A, > 0. Here, we only consider 
three dimensionally expanding genuine voids, > Ai > A2 > A3 and 5 = X^j=x ^^ ^ ^■ 

Now, under the premise of the strong correlation between the void shape, Ijj, and the 
underlying tidal tensor, Tij = q q , we have the following relations between the semimajor 
axis length, pj, and the eigenvalues of the tidal tensor, A^: 
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Figure 2. Wc show the ellipticity PDF for Rl = 5 and 10/i~^Mpc at z = 0. PDF is normalized so 
that the peak values are 1. Solid lines show T = 0.5 and stars show the peaks of each PDF. 



where we assume Ai > A2 > A3. Using these, one can define the PDF of e and r] as 



where /x = 1 — 77 and z^ = 1 — e, and 
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From eq. (2.4), we see that the probability of having two of the three eigenvalues, Aj, 
being equal is suppressed due to the last three factors (i.e., (Ai — A2)(A2 — A3)(Ai — A3)). 
Therefore, we have small probabilities of having ellipsoids with special configurations such 
as e = ?7 (prolate spheroids), rj = (oblate spheroids) and e = (sphere): these special 
configurations correspond to Ai = A2, A2 = A3 and Ai = A3, respectively. In figure 2, we plot 
the elhpticity PDF of eq. (2.11) for Rl = 5 (left) and 10 h'^Mpc (right), where <5^, = -0.9. 
PDF is normalized so that the peak value is unity (stars). Note that for a fixed density 
contrast, 5, the exponential factor in eq. (2.4) is maximized when Ai = A2 = A3 allowing 
nearly spherical voids to have a finite probability. 

We also plot the line separating the prolate and oblate ellipsoids (i.e., T = 0.5) (solid 
line). We see that locations of the ellipticity PDF peaks lie roughly on the line of T = 0.5, 
and PDF shifts toward smaller ellipticity both in e and rj for larger voids (i.e., smaller cr/j^ 
for a fixed set of cosmological parameters). 
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Figure 3. We show the 232 identified voids from the SDSS DR5 based on the catalog of Foster 
& Nelson [8], and the corresponding ellipticity PDF. From the catalog, we find e = 0.33, fj = 0.21, 
T = 0.68, Sy = —0.97, z = 0.11 and Rl = 7.3/i~^Mpc. PDF is normalized so that the peak values 
are 1, and the solid line shows T = 0.5. Note that the PDF is drawn only as a reference, and we do 
not correct for bias and redshift space distortion. 



Also note that the ellipticity PDF in the real space has almost no preference for prolate 
or oblate ellipsoids (T ~ 0.5). For Rl = 5/iMpc~ , 6 = —0.9 and z = 0, we have mean 
ellipticities, e and f] as follows, 
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and for Rl = l/iMpc~\ 6 = -0.9 and z = 0, we have e = 0.62, r? = 0.27 and T = 0.55. 

However, previous observations and numerical studies indicate a preference for prolate 
ellipsoids. From voids detected in a ACDM N-body simulation. Platen et al. [31] claimed 
the axis ratio of best-fit ellipsoids to be 1 : 0.7 : 0.5 (i.e., e = 0.5, r] = 0.3 and T = 0.68), 
and FN09 also found a preference for the prolate ellipsoids from the voids detected in the 
SDSS DR5 with e = 0.33, r] = 0.21 and T = 0.66. In figure 3, we show the distribution of 
ellipticities of 232 voids from the SDSS DR5 based on the catalog of Foster &: Nelson [8]. We 
also show the corresponding ellipticity PDF just as a reference. The volume limited sample 
has a total of 52281 galaxies in a volume of 0.021/i~^Gpc^ {fig ~ 2.45 x 10~^/iMpc~^). In the 
e-i] plane, we clearly see the dominance of prolate voids, while analytic PDF has a slight to no 
preference of prolateness in the shape of voids. We think that the difference arises from the 
way we calculate ellipticities of voids. Earlier works extracted void shapes by calculating the 
shape tensor from volume elements inside detected voids. In this way, an extracted ellipticity 



is sensitive to the boundary shape of the void as defined from a given void finder. Since 
a boundary is defined from galaxies in an overdense region, 5 > 0, the tidal field around 
the boundary is subject to a stronger distortion from the external gravitational force. An 
unconditional ellipticity PDF of Doroshkevich [7], on the other hand, describes a tidal field 
at and around the density minima. Therefore, when we compare the analytic ellipticity PDF 
of PL07 against that of a simulation or galaxy survey, it is important to know whether the 
extracted ellipticity truly represents the shape of the tidal field around density minima. 

3 Linear Galaxy Bias in the Real-Space Void ellipticity PDF 

In the previous section, we have related the ellipticity PDF of voids to the tidal field generated 
by the matter density distributions. However, in reality, we obtain a void ellipticity PDF 
from a given galaxy distribution. In order to deduce the underlying matter density field, 6m, 
from galaxy distributions, we use the well known linear galaxy bias, 6i, such that 5^ = hL&m- 
Here, 6y represents the density contrast of a void calculated from the number density of 
galaxies within the void, nyg, and that of the cosmic mean, n^, (i.e, average number density 
of galaxies in a survey volume) , 

6, = !!Z£^^. (3.1) 
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Since the unconditional PDF of Doroshkevich [7] gives eigenvalues of tidal tensor, Tjj, gen- 
erated by the matter density field, we use 5m = ^v/^l in Eqs.(2.8) ~ (2.13), where 6^ is a 
linear galaxy bias to be calibrated from an observed ellipticity PDF. The filtering scale, R^, 
should also be calculated following the definition of the linear bias, 

R^ = R^{l + 5,/bLfl\ (3.2) 

Here, an Eulerian size of void traced by the matter density field, i?^, is not a direct observable 
in a galaxy survey (i.e., the void size measured by the galaxy distribution, i?^, is a direct 
observable). Patiri et al. [28] showed that the radial density profiles of voids defined from 
dark matter and halo distribution of N-body simulations have similar shapes, and at large 
radii, haloes trace the dark matter. Based on this result, we set R^ ~ i?|; here, and have 

R^ = Rl[l + 5JhLfl\ (3.3) 

Figure 4 shows the effect of linear galaxy bias on the real space ellipticity PDF both in 
1-D (bottom right) and in 2-D (top left to bottom left). We normalize the 2-D PDFs to their 
peak values, and set the parameters to Rl = 5/i~^Mpc, 6y = —0.9 and z = 0. In the 1-D 
plot, we have PDF for four different linear galaxy biases, b^ = 1.0 (solid), l.l(dotted), 2.0 
(dashed) and 5.0 (dot-dashed), and in the 2-D plots, we have 6^ = 1.0 (top left), 1.1 (top 
right) and 2.0 (bottom left). 

In the 2-D plots, the solid lines show T = 0.5, and the stars show the peaks of each 
PDF. Regardless of the values of bias, peaks of PDF lie on the solid lines (T = 0.5), and the 
ratio of prolate to oblate spheroids is ~ 0.5 as in the case of unbiased (i.e., b^ = 1) ellipticity 
PDF. 

We clearly see an effect of bias on the void PDF as we increase the galaxy bias bi from 
its unbiased value, 6l = 1, to biased value, bi > 1. As we increase the bias from 6^ = 1, 
the ellipticity PDF peaks at smaller ellipticities more sharply. For bi > 2.0, the shape of 
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Figure 4. In this figure, we sliow tlie effect of linear galaxy bias on the real space void ellipticity 
PDF both in 1-D {bottom right) and 2-D {top left to bottom left). The 2-D PDFs are normalized so 
that the peak values are 1. In the 1-D plot, we have PDF for four different biases, 6l = 1.0 (solid), 
l.l(dotted), 2.0 (dashed) and 5.0 (dot-dashed). In the 2-D plots, the solid lines show T = 0.5, and 
the stars show the peaks of each PDF. Here, we have used Rl = 5h^^Mpc, Sy = —0.9 and z = 0. 



PDF converges quickly, and becomes insensitive to the further change in b^. As a larger 
bias reduces the matter density contrast, 5m, for a fixed value of 5^, a smaller matter density 
contrast yields a less pronounced tidal field around the density minima, yielding smaller 
ellipticity. 

In Eqs.(2.8) ~ (2.13), we see the effect of bias through a^^. Since the Lagrangean size of 
void, R^, (as defined in eq. (3.3)) sets the smoothing scale of the density field for calculating 
cjij^, the larger the K^ is, the smaller the 0"^^ becomes. Again, a smaller rms fluctuation 
of the matter density field leads to a weaker tidal field, and hence, the smaller ellipticity of 
the void. We see that, although the linear bias, bL, and ag, which is directly proportional 
to the ajij^ , are nearly degenerate in the parameter space of the biased ellipticity PDF, their 
correlation is directly opposite in the case of the galaxy power spectrum, Pg{k) oc b\a'^. This 
feature of the biased void ellipticity PDF allows us to better constrain both the linear bias, 
bL, and the rms density fluctuations, a%, when combined with the galaxy power spectrum. 



4 Void in the redshift space 

In order to compare the analytic and observed ellipticity PDFs, we need to understand the 
effect of the redshift space distortion on the elhpticity of a given void. UnUke galaxies and 
galaxy clusters, which reside in density maxima, voids reside in density minima, where the 
diverging velocity field pushes the boundary and the field galaxies inside the voids outward, 
expanding their volume in all directions [14]. 

A real-space inertia tensor is defined in eq. (2.1), and that of the redshift-space is 
similarly defined as 

/*■ = 'S^ r" r" ■ (A l\ 

where A^vg is the number of void galaxies, and x^ is the redshift-space position of an a-th 
void galaxy. An inertia tensor is symmetric by definition (i.e., lij = Iji), and traces a shape 
of a void from an underlying void galaxy distribution. Park &: Lee [27] identified voids from 
the Millennium simulation using the void finding algorithm of Hoyle &: Vogeley [12], and 
extracted the ellipticity PDF of voids using the inertia tensor using void galaxies as tracers. 
They found that the distribution of ellipticities of voids from the Millennium simulation 
follows that of linear theory (Zel'dovich Approximation) prediction in real space. 

Here, we study the effect of the redshift space distortion on the ellipticity PDF. First, 
we calculate the redshift space inertia tensor as follows. 

A mapping between the real-space and the redshift space is given by 

(4.2) 
(4.3) 
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where / = ^\^^ is the growth rate, and the unit vector, Xg, is along the line of sight. Here, 
u = m, where v is the peculiar velocity of a void galaxy, and T-L{a) = aH{a) is a comoving 
Hubble rate. 

We define k(x^) to quantify the strength of the redshift space distortion as 

.(x^) = -M:^. (4.5) 
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With k(x'"), we re- write eq. (4.4) as 

xl = xl{l + k{^^)\ (4.6) 
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and we relate the redshift-space inertia tensor to the real-space inertia tensor as follows, 
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where 



J fJ-aUa 



and /i is the angle between the peculiar velocity, u, and the line of sight (i.e., ^ = u • £3). 

Once we have the redshift space inertia tensor either numerically, from eq. (4.1), or 
analytically in terms of k, from eqs. (4.7) ~ (4.12), we can calculate the ellipticity of a given 
void, e and rj, from eigenvalues of the redshift space inertia tensor, (If > /I > I§) as in eq. 
(2.2). 

Since k is a function of the peculiar velocity field, we can use linear theory to calculate 
K for an arbitrary density field, S^x^), as, 

/xVu(xO _ /xVV<I>(x'-) 



"3 



AnGa^px^ 



^ ^ '^'^rik-iS|.(k)e*^'''' 



iTTGa^pxlJ (27r)3 

xlJ (27r)3'^ ""^ k^ " 

_f_ f dk .^ f ^^^ f ^^^^.^^i^vikx-- 
4io (27r) 



f°° dk /■^'^ f^ 

Jo [^'^r Jo J-i 



(4.14) 



where fik = cos^^. Here, we used the linear relation between a peculiar velocity and a 
gravitational field derived from Zel'dovich Approximation [13], 

^■^ V^, (4.15) 



AirGa'^p 
and a Poisson equation, 

V^^ = A-KGa^p6. (4.16) 

We explicitly write a cosine angle between k and x, in a spherical coordinate as 

cos 7 = cos 6k cos9x + cos{(j)k — <i>x) sinSfc sinS^;, (4-17) 
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where k • x = kxcos'j. Here, in order to calculate the effect of peculiar velocity on the 
void shape, we need to know the density profile of a void, 5(x) in the real space. Sheth &: 
van de Weygaert [36] studied an evolution of a radial density profile of an isolated void, and 
found that the density profile evolves into an increasingly similar shape to a top-hat function. 
With N-body simulations both in sCDM and ACDM, Colberg et al. [3] found that the radial 
density profile of a void is universal, having p{r < res) / p{refi) oc exp[(r/refr)°], where rgfj 
is the effective radius of the void and q ~ 2. These analytic radial density profiles are also 
seen in the voids identified from the galaxy distribution of the Sloan Digital Sky Survey as 
a sharp rise of density contrast near the boundary of the void [26, 28]. In this section, we 
approximate the radial density profile of the void, d{r), to be a top-hat function for simplicity, 

^, r\ _ j ^ within the void , . 

^^""^-{o outside the void ' ^ ^' 

where 6 is the average density contrast within the void. We will discuss further about more 
general density profiles at the end of this section. 

Below, we will show analytic solutions for the redshift space ellipticity for some of the 
limiting cases. 

4.1 Spherical Void 

First, we consider a sphere with zero real space ellipticity, e'' = and rj^ = 0, with three 
eigenvalues of the inertia tensor of I[ = I2 = -^3- For a spherical top hat density field, we 
have 



6(k) = AttWS-^ — '-. (4.19) 

kR 



Therefore, eq. (4.14) becomes 



_ 47rfi3/J p dk . JijkR) 
''^'' ^ " 4 io (2^)' kR 

(•2-K rl 



X I d(t>kj dpkf^ke"'-'' . (4.20) 

We first integrate the angular portion of the integration, f dilfc/i^e , by expanding the 
exponent into Legendre polynomials, 



gik-x- 



J^(-i)'(2/ + l)ji{-kx)Pi{cosj), (4.21) 

I 

and using the Spherical harmonics addition theorem (a.k.a., Legendre addition theorem). We 
find 

/ dnk^^k(^^'''' = 47riji (kx) cos e^. (4.22) 

Finally, we obtain k for a spherical void as 

K(x^) = — ^y^ ^^^k^^[4.^nikx)cos0.] 



2R^f6cose.^ ^ TTx \ _ f6 
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As we see from eq. (4.23), k. is independent of positions in the void, and therefore the three 
eigenvalues of the redshift space inertia tensor for a redshift space void, /f — -^1 — -^1' ^^^ 

If = /[(I + k)2 (4.24) 

II = II (4.25) 

/| = II, (4.26) 



and the ehipticity, e and rj are 



1 



1 + K 1-/5/3' 



(4.27) 



foUowing eq. (2.2). Now, the hnear growth rate, / = ^^^ ^^' , for the ACDM model is 
approximately given as [11, 29], 

f{z) = n{zfl\ (4.28) 

where 

For Vtrn = 0.25 and Oa = 0.75, / ~ 0.5 at z = 0, and thus, for a spherical void with an average 
density contrast of (5 = —0.9, we have e* = r/* ~ 0.13. Since there is no special orientation for 
a spherical void, any void with sufficiently small real space ellipticity, e*" ~ r/^ ~ 0, becomes 
a prolate ellipsoid (i.e., e* ~ ry'^ 7^ 0) in the redshift space with its ellipticity given by eq. 
(4.27). 

4.2 Spheroidal Void 

When two of the three lengths of the semiaxis are equal, we have a spheroidal void. We set 
the direction of the longest and shortest semiaxis to be along the line of sight for prolate 
and oblate spheroids respectively. As in the spherical void, we assume a spheroidal top hat 
density contrast and solve eq. (4.14) to have k to the first order in e*", 

K(x^•)c^-^(lT2e^) (4.30) 

where "— " and "+" signs are for prolate and oblate spheroids respectively. Here, we set the 
longest axis along the line of sight for prolate spheroids, and set the shortest axis along the 
line of sight for oblate spheroids. We compared the analytic prediction against numerical 
calculations, and found that 

«:(x^):^-^(lTe^) (4.31) 

fits the result better. For each case, we have redshift space ellipticity as a function of real 
space ellipticity and k as follows, 

.• = 1 - i^. (4^32) 

e^ = l-(l-e^-)(l + K). (4.33) 
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Figure 5. We show the spatial distribution of the value of k(x'") in x-z plane for prolate spheroids 
with its longest axis along the line of sight (top), and oblate ellipsoids with its shortest axis along the 
line of sight {bottom) for different real space ellipticitics, e''. We also plot real space (solid lines) and 
rcdshift space (dotted lines) shapes of voids. Here, we use / = 0.5 and S = —0.9. 



In contrast to the spherical void, whose redshift space ellipticity is always larger than that 
of real space, ellipticity of a spheroid either increases or decreases in the redshift space. 

In figure 5, we show the result of numerical calculations for the spatial distribution of 
the value of fi;(x'^') in the x-z plane for prolate spheroids with its longest axis along the line 
of sight (top), and oblate ellipsoids with its shortest axis along the line of sight (bottom) for 
different real space ellipticities, e''. In the figure, the center of the spheroids corresponds to 
the origin of the coordinate, and the x-z plane cuts the middle of the spheroid. We also plot 
the real space (solid lines) and redshift space (dotted lines) shapes of voids. Here, we use 
/ = 0.5 and 6 = -0.9. 

We see that the deformation of the void shape, k ex U3, is the largest when the shortest 
semiaxis is along the line of sight due to the largest line-of-sight peculiar velocity. This result 
is in agreement with the earlier work of Icke [14]. When the longest semiaxis is along the 
line of sight (i.e., top figures), ellipticity only increases in the redshift space, while ellipticity 
only decreases when the shortest semiaxis is along the line of sight (i.e., bottom figures). For 
the spherical void, e^ = 0, we have e** = 0.13 as expected from eq. (4.27). For prolate and 
oblate spheroids with e^ = 0.2, 0.4 and 0.6, we have k = 0.12, 0.09 and 0.06, and k = 0.18, 
0.21 and 0.26 from the numerical calculations respectively. This is in a good agreement with 
the fitting function of eq. (4.31): for prolate and oblate spheroids with e'' = 0.2, 0.4 and 0.6, 
eq. (4.31) yields k = 0.12, 0.09 and 0.06, and k = 0.18, 0.21 and 0.24, respectively. 

Using Eqs.(4.33) and (4.33), for prolate and oblate spheroids with e^ = 0.2, 0.4 and 0.6, 
we have redshift space ellipticities of e'^ = 0.29, 0.45 and 0.62, and e* = 0.06, 0.27 and 0.50 
respectively. In general, the fractional change in the ellipticity, (e* — e^)/e^ , becomes smaller 
for a larger e''. 

Although we have constant n here (i.e., independent of the position in the void), this 
is not true for arbitrarily oriented voids. In the case of tilted voids, k;(x'^') is a function 
of the real space coordinate, x*", and has both positive and negative values by definition, 

k(X) CX U3/X3. 
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Figure 6. Wc show the effect of redshift space distortion on real space void cUipticities, e^ and ry'". For 
a given set of ellipticities in real space, we plot the change in ellipticities, Ae(e'', ry*") = e'*(e'", r]^) ~ e^ 
{left) and Ar]{e^,7f) = jf{e^,ri^) — rj^ (right). Here, we have assumed a top-hat density profile with 
d = -0.9 and / = 0.5. 



4.3 Redshift-Space EUipticity Distribution Function 

For a given ellipticity in real space, the corresponding redshift-space ellipticity can either 
increase or decrease, depending on the real space shape and orientation of the void with 
respect to the line of sight, except for a spherical void, whose ellipticity can only increase. 
Here, we calculate changes in ellipticities for a given set of real space ellipticities, e^ and r/^, 



Ae(e^7?'■)EEe^(e^r?^)-e'• 



(4.34) 
(4.35) 



averaged over possible orientations. We numerically solve eq. (4.14) for k(x'"), with a top- 
hat density contrast, 5(x'"), and calculate the redshift-space ellipticity, e'* and r/**, from the 
redshift-space inertial tensor. If-. In figure 6, we show Ae(e^,r7^) {left) and Arj{€'^,ri^) (right) 
for 5 = —0.9 and / = 0.5. As expected, a spherical void in the real space becomes a prolate 
spheroid (i.e., e = 77 / 0) in the redshift space with its ellipticity given by eq. (4.27). Also, 
we see a clear preference for positive Ae for any triaxial void, especially at a small e^. As 
for Ar], spheroidal voids tend to be triaxial voids in the redshift space. From both Ae and 
Ar], we find that an oblate void (T < 0.5) in real space tends to be a more prolate shape 
(r > 0.5) in the redshift space, and we have more prolate voids in redshift space ellipticity 
PDF. 

In figure 7, we compare the ellipticity PDF of real space and redshift space for Rl = 
4h~^Mpc (top) and Rl = 9/i~-'^Mpc (bottom). Here, we use hL = 1 and 5 = —0.9, while 
varying the linear growth rate, /, from 0.5 to 1.5. We show the redshift space ellipticity 
PDF together with that of real space as a reference. Since the expansion rate of the void is 



15 



10.0 



Q- 1.0 r 



0.1 
10.0 



Q- 1.0 r 



0.1 



0.0 



^^.xs,^ R =4 Mpc/h 




^ / ■ 







R =9 Mpc/h 



» \ real space 
f = 0.5 
f=1.0 
f=1.5 



0.2 



0.4 0.6 

ellipticity (e) 



0.8 



1.0 



Figure 7. We compare the ellipticity PDF of real space and redshift space for Rl = 4/i^^Mpc {top) 
and Rl = 9/i^^Mpc {bottom). Here, we use bL = 1 and S = —0.9, while varying the linear growth 
rate: / = 0.5 (dotted), 1.0 (dashed) and 1.5 (dot-dashed). We also show the real space ellipticity 
PDF (solid) as a reference. 



set by the matter density distribution inside the void, in linear theory, it is proportional to 
the linear growth rate. As the effect of the redshift space distortion on the ellipticity PDF is 
more prominent for a larger f{z) and smaller e^' , the redshift space ellipticity PDF becomes 
more skewed in shape with the mean ellipticity being higher than that of real space, e* > e^ . 
Also note that the left tails of each redshift space ellipticity PDF approximately show the 
redshift space ellipticity, e'*, of a spherical void, e^ = 0, and thus, we have 



^left ^s { r 

e ~ e (e 



0) 



1 



1 



3 



(4.36) 



Therefore, in principle, with only the location of the left tail of PDF determined, we can put 
a constraint on f5. This is a great benefit of using morphological information of structure, 
which is directly distorted in the redshift space. Although what we can actually measure 
is f6 = j35v, where j3 = f /hi, the entire shape of ellipticity PDF can lift the degeneracy 
between / and hi. We can easily see this in the figure 4, where the increase in bi decreases 
the mean ellipticity, e^, shifting the right tails of PDF toward lower e*", while the left tails are 
almost unaffected. Unlike a galaxy power spectrum, P{k), which gives a constraint on the 
linear growth rate that is degenerate with the linear bias, /3 = f/bi, a void ellipticity PDF 
has a potential to constrain the linear growth rate separately from the other cosmological 
parameters. In other words, for a given galaxy survey, we have the situation in which we 
"6ii?/ one, get one free" . 
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4.4 Void With General Radial Density Profile 

So far, we have derived the redshift space ehipticity PDF for the top-hat density contrast, 
(5(x); however, the reaUstic radial profile of void density contrast does not necessarily follow 
the top-hat shape. The main parameter controlling the shape of the redshift space ellipticity 
PDF is k(x'^') as in eq. (4.14), and its value is given by k(x'^') = A{x'')f6, where ^(x'') is 
solely determined by the density profile of a void (e.g., a void with a spherical top hat density 
profile has ^(x'') = —1/3). 

As an example, we parametrize the radial density profile of the void following the 
accumulated density profile of voids from the N-body simulation of Colberg et al. [3], p{< 
r)/p{rcs) oc exp[(r/rcff)°], where the dentity profile at a given radius, p{r), is given as follows 



p{r) ex exp 



r \ 
rcsj 



a I r 



3 yvcs 



(4.37) 



Here, instead of using the original form of Colberg et al. [3], we use an exponential 
density profile, 

6{r) = 6oi ^zr^l 1' (4.38) 

where we forced p{rcf[) = P, at the effective radius, reS, (i.e., 6{rcs) = 0), and 5q = 6{r = 0). 
Figure 8 shows the profiles of the density contrast, S{r), for top-hat (solid) and exponential 
profiles of eq. (4.38) with a = 1 (dotted), 2 (dashed) and 4 (dot-dashed). Note that for a 
top-hat density contrast, we have 5 = 6o, while for the exponential profiles, we have 



/o'^« r^6{r)dr 
L °^ r'^dr 



(4.39) 



where 6{r) is given as eq. (4.38). In figure 9, we show the redshift space ellipticity PDF for 
/ = 0.5 (green), 1.0 (red) and 1.5 (orange) together with the real space ellipticity PDF (blue). 
Here, we use Rl = 9/i~^Mpc and the radial density profile of the density contrast is set to 
be either a top-hat (solid), 6{r < Rmid) = S, or an exponential (dotted) given by eq. (4.38) 
with Teff = -Rvoid and a = 2. We use 6o = —0.40 and —0.82 for the top-hat density profile 
and the exponential profile, respectively, in order to keep the same mean density contrasts, 
5 = —0.4, for different density profiles. Here, we clearly see degeneracy between the shape of 
the potential and the linear growth rate, /. Therefore, a better understanding of the density 
profile of a void is required when constraining the linear growth rate from the redshift space 
ellipticity PDF. 

5 Effect of Void Finding and Poisson Noise on the Void Ellipticity PDF 

In the previous sections, we have shown how to incorporate the galaxy bias and redshift space 
distortion into the calculation of the ellipticity PDF. In this section, we investigate the biases 
and errors that arise when using the inertia tensor, /jj, as an estimate of the void ellipticity. 
Since a void is a sparse region in a galaxy distribution, the number of galaxies inside a void, 
A'vg, is limited. Therefore the error in the ellipticity calculated from the inertia tensor of the 
void galaxy distribution is dominated by the Poisson statistics. 
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Figure 8. We show profiles of the density contrast, 5(r), for top-hat (solid) and exponential profiles 
of eq. (4.38) with a = 1 (dotted), 2 (dashed) and 4 (dot-dashed). Here, we assumed that the density 
at the effective radius, roff , becomes that of the cosmic mean, p(rcs) = p, or, ^ivcs) = 0. 
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Figure 9. We show the redshift space ellipticity PDF for / ~ 0.5 (green), 1.0 (red) and 1.5 (orange) 
together with real space ellipticity PDF (blue). Here, we use Rl = 9/i~^Mpc. The radial density 
profile of density contrast is set to be either a top-hat (solid), 6{r < i?void) = S, or an exponential 
(dotted) as in eq. (4.38), with VeS = Rmid and a = 2. We use 6o = —0.40 and —0.82 for the top-hat 
and the exponential density profile, respectively, in order to keep the same mean density contrasts, 
S = —0.4, for different density profiles. Here, we see degeneracy between the shape of the potential 
and the linear growth rate, /. 



Let us define two stochastic variables, x and y, as follows, 

X = Hx + Ax, 
y = l^y + Ay, 



(5.1) 
(5.2) 



where Hx and i-iy are the means of x and y, respectively, and (A.x) = (Ay) = 0. We have 
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variances and covariance defined as, 



al ^ {{X - ^,.f) = {{Axf) , 

(Txy = {{X - /ix)(y - /-fy)) = {AxAy) . 



A mean of a rational is given by the following form, 

_ / fix + Ax 
^''/y - \ M, + Ay 

Expanding Ay up to the second order, we have 

Ax 



l^x/y 







a. 



xy 



f^y f^xfJ-y 



Similarly, variance of the rational is 



(5.3) 
(5.4) 
(5.5) 



(5.6) 



(5.7) 



a: 



x/y 



' fix + Ax ^ 
jiy + Ay^ 

^i^i^ + 'i 

" ,2 ,,2 ^ ,,2 



l-l^x/y 



2a. 



xy 



f^y \ f^x f^y fJ-xfJ-y 



(5.8) 



Now, ellipticity, e, is defined ase = l — — = 1 — \/ -r-, where pi > P2 ^ Ps are the three 



pi 



principal axes of the ellipsoid, and h > h > h are the eigenvalues of the inertia tensor lij. 
From pf (X li, we have 



Api _ l AIj 



and therefore, we have 



and 



l^e = l- f^p,/p, ^ 1 






1 + 



O", 



pi 



cr. 



Pip-i 



f^Pl /^P3 






(5.9) 



(5.10) 



^2^ 2 ^^^llk + 

"^ "Pz/Pi - 1,2 I ,,2 ^ 



^P3 



^<7piP3 



l^K + 4 



fJ'pi l^pz I 



(5.11) 
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From eq. (2.1), we have 






1 

- y^Xio^Xj,^, yO.iZ) 

and in the continuous density field, we have 

lij = d X n{:x.)xiXj 

(fx [n + 5n(x)] XiXj = fij.. + A/jj, (5.13) 



/ 



*void 

where 

^Li^^ = I (fx nxiXj, (5-14) 

Alij = / d^x 6n{:x.)xiXj. (5.15) 

Also, the variance of the inertia tensor is given as 

<, = ((A/.,)^> = (4>-4r (5-16) 

5.1 Void With a Well-defined Boundary: {pi) — ;• Pi 

Let us calculate the bias in the mean and the variance of the ellipticity calculated from a 
rectangular void with a diagonalized inertia tensor of Ji = In, I2 = I22 and I^ = 133. The 
number of void galaxies is N^g. Here, we assume an idealized setup, where the boundary 
shape of a void is well defined (i.e., [pi) — )• Pi) to see whether the estimator of an output 
ellipticity based on the inertia tensor yields unbiased result. For more realistic estimate on 
the bias in mean and the variance of output ellipticities, see § 5.2. 

For a rectangular void with volume, Vyoid = 8PiP2-f3 (i-e., Xi G [—Pi-, Pi] V i), where Pi 
is a well-defined length of the principal axis of the rectangular void, we have 

, = I cfix nx\ = ^nV.oidPf = \N.^Pf , (5.17) 



*void 



and 



al = / / d^xd'x' (n(x)n(x')) x^xf - ^ 

JV^oid JV ,7 

voza ' void 

= / d^xd^x' n5D{y^-y^)x'ix'i - ix]. 



= n d^x xf - nj_ = -nV^oidPt - A 

V ^ 5 

= ^^vg^' - ^^vg^'> (5.18) 

where we have assumed that the distribution of void galaxies follows the Poisson statistics 
with (n(x)) = n and (n(x)n(x')) = n6Di^ — x'). Similarly, we have 

crijj = n d^x xjx] - fii^fii^ 



l-N,,P^Pf - \NlPfPl (5.19) 
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Therefore, Eqs(5.10) and (5.11) are re-written as follows, 

^- = f 5^- (^■^'> 

We see that for a given rectangular void with N^oid = 20, we can accurately extract its 
ellipticity from void galaxies to ~ 1%, but its standard deviation remains significant (i.e., 
a^ ~ 0.14 for a void with zero ellipticity). As is expected, the variance of ellipticity, a'^, scales 
as inversely proportional to the number of void galaxies, iVvg- Note that the above result is 
derived in the limit of {pi) — )■ Pi. 

5.2 Void With an Undefined Boundary: (pi) / Pi 

With a small number of field galaxies randomly distributed inside the well-defined boundary, 
Pi, we generally have {pi) ^ Pi. For example, let us imagine a spherical boundary, where 
Pi = P2 = P3 and ejn = 1 — -^ = 0, with N^g = 4 field galaxies. This spherical boundary 
represents the shape of the underlying tidal field, and we sample this by 4 galaxies. No 
matter how we distribute those four galaxies, unless the four galaxies align at each of four 
apexes of an equilateral tetrahedron, calibrated ellipticity cannot be zero, eout / e^n- As a 
result, we have a significant deviation of (cout) from its boundary ellipticity, ej„ = 1 — -^, 
and a correspondingly large standard deviation. 

As we increase A'vg, those randomly scattered particles fill in the bounded space more 
densely, and (eout) gradually approaches the boundary ellipticity, e^- This uncertainty in 
€out due to small A^vg introduces an additional bias and scatter to the ones given by eqs. 
(5.20) and (5.21). 

For the case of elliptical voids, the analytic prediction of the mean and variance is 
non-trivial; thus, we have performed numerical experiments. First, we create an ellipsoidal 
boundary with a fixed size and ellipticity (i.e.. Pi, P2 and P3), and then fill this bounded 
volume by N^g random particles. For each realization, we calculate the inertia tensor and its 
eigenvalues (pi, p2 and ps), where the output ellipticity is defined as 



. -1 ^3 
Pi 




(5.22) 


Jr -L 

Vout = 1 , 

Pi 




(5.23) 


while the input ellipticity is defined as 






e -1 ^' 
'*"- Pi' 




(5.24) 


1 ^2 

r?m = 1 p . 




(5.25) 


Finally, we have a bias in the mean and variance of the output ellipticity, e°"* 


as follows 


/ \ ^3 - 
\^out) ^in — J-, \ 
Pi 


/P3\ 
\Pl/' 


(5.26) 


^Lt = (^L) - < 


,^out) ■ 


(5.27) 
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Figure 10. Bias and variance in the measured ellipticities. We place A^vg particles within a boundary 
of a given set of size and input ellipticities, e™ and ry™, and calculate the mean ellipticities, (e°"') 
and (t?°^*), and their standard deviations a^ and ar^, from 500 realizations for each N^g and input 
ellipticity. Top figures show the biases in finding the input ellipticities (i.e., true underlying ellipticities 
of a tidal field), (e°"*) — e™ and (77°"*) — 77™, and bottom figures show the standard deviations, tJe and 



Figure 10 shows the result of 500 realizations for a given set of e"^ and -/Vvg- The 
top figures show the bias in the mean of the ellipticities e (left) and r] (right), and the 
bottom figures show the standard deviations of e (left) and r] (right). We see that both bias 
in the mean and variance of eout are large enough to change the shape of ellipticity PDF 
significantly for A^vg ^ 100. In order to correct the ellipticity PDF, we define a response 
function, i?(eout|eiin -^vg), such that 



p(eout|o-i?,i,,<J,A^vg) 

= / deinR(eonthn,N^g)p(em\(^RL,S)- 



(5.28) 



In figure 11, we show the normalized histograms for the number of voids with output ellip- 
ticity, Eout; from the numerical simulation with a given A'vg and ein- We also show Gaussian 
distributions (red lines) with the mean (eout) and variance cf'^^^^, 



R(^out\^m, N, 



1 



vgy 



V2^ 



: exp 



(Cout — (Cout))^ 



(5.29) 



where both mean and variance are functions of the input ellipticity and the number of void 
galaxies (i.e., (eout) = (eout(ein, A^vg)) and a^^^^ = a^^^^(ein,N^g))- In figure 12, we show the 
ellipticity PDF without the correction for the response function (black lines) as well as the 
corrected ellipticity PDF (red lines) for four different redshifts, z = (solid), 1 (dotted), 2 
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Figure 11. We show the normahzed histograms of output ellipticities, eout for a fixed input ellipticity, 
Ein, and a number of void galaxies, N^g. Wc find that the Gaussian profile of eq. (5.29) with mean, 
(cout), and variance, cTeo^t, fits the measured histograms (red lines) well. Also, we show the locations 
of the input ellipticity, £;„, with vertical lines. 



(dashed) and 3 (dot-dashed) with different values of A^vg- Here, we use Rl = 4Mpc~^ and 
^m = —0.9 to calculate the original PDFs at four different redshifts. We see that the shapes 
of PDF with higher redshifts are more susceptible to the small iVyg effect. For a typical 
void studied in PLOT using the Millennium Run catalog at z = 0, N^g is roughly between 
20 to 60, Rl ranges from 4 to 6.4/i~^Mpc, and 5^ ~ —0.9. We clearly see that voids with 
small initial ellipticity ejn are most vulnerable to the Poisson noise, and with A^vg ~ 10, the 
original shape of ellipticity PDF is almost washed out to leave only the distribution of the 
Poisson noise. On the other hand, the ellipticity PDF at low redshift {z = 0) maintains its 
original shape with the number of void galaxies as small as A'vg ^ 45. In short, a shape 
with small ellipticity is hard to trace with a small number of tracers. For example, if we 
have a sphere of radius R = 10/i~^Mpc in the simulation box with a resolution of l/i~^Mpc, 
the possible number of configurations of randomly scattered tracing particles with A^vg = 4 
will be 4187C4 ~ 1.3 X 10^^, while the possible number of configurations of an equilateral 
tetrahedron (i.e., e = 0) is vanishingly small. Also, the scatter of the output ellipticity a^ 
will not center around ein for a small ellipticity as < e < 1, and as a result, (eout) converges 
to ein only slowly as we increase A'vg- On the other hand, the ellipticity of a highly elongated 
shape can be well traced even with a limited number of tracers. 

For A^vg ^ 100, we use a fitting function derived from figure 10, 






= ai(ei„)Af-/i(^-) 



(5.30) 
(5.31) 



where ai, a2, /3i, and /32 are given as polynomials to the third order, Po+pi^m+P2^in~^P3^tn^ 
and we show the best-fit coefficients, po, pi, p2 and ps in table 1. 
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Figure 12. Wc show the real space eUipticity PDFs without Poisson noise (black lines) at z = 
(solid) , 1 (dotted) and 2 (dashed) . Wc then convolve each noise- free PDF with the response function 
given by eq. (5.29), and show the resulting PDFs for four different A^vg (red lines). 





Po 


Pi 


P2 


P3 


ai 


1.282 


2.384 


-7.280 


3.587 


02 


0.2636 


-0.3103 


0.9557 


-0.9404 


Hi 


0.4500 


1.708 


-0.8460 


-0.02668 


P2 


0.3650 


-0.6260 


1.680 


-0.9296 



Table 1. We show the best-fit coefficients, po, pi, p2 and p^ for ai, a2, /3i, and /32 of eqs. (5.30) and 
(5.31), where ai, a2, /3i, and (32 are given as polynomials to the third order, po +Piein +P2^fn +P3^fn- 



6 Voids from N-body Simulation 

To test the analytic predictions made in the previous sections, we use a catalog of galaxies 
from the Millennium simulation [5, 6]. In a comoving box of Lbox = 500 /i^^Mpc, there are 
26690265, 26359329, and 23885840 galaxies in the simulation at z = 0, 1, and 2, respectively. 
For each redshift, we select a subset of galaxies whose halo masses are above Mh = 2.6x10^^, 
2.6 X 10^2^ and 8.6 x W^^h'^MQ. We find 182081, 121454, and 261471 such halos for z = 0, 
1, and 2, respectively. 

Here, we use the void finding algorithm of FN09 to identify voids from the subsets of 
galaxies. Their algorithm first divides a given set of galaxies into wall galaxies, which define 
boundaries of holes, and void galaxies, which are allowed to be within holes by calculating 
the distances to the third nearest galaxies, D3. When a galaxy has the third nearest distance, 
I?3, smaller than the user-defined threshold distance, Rs > D^, the galaxy is labeled as a wall 
galaxy, and when R3 < D3, the galaxy is labeled as a void galaxy. Here, R3 = {D3) + Xctd.^, 
where {D3) and a^a are the mean third distance and its standard deviation, respectively. 
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and A is a user-defined parameter. We use the value of A = 2 recommended by FN09. Then, 
a hole is defined as a sphere with maximal radius in a galaxy distribution containing no wall 
galaxy, and by definition, a hole can contain void galaxies. Finally, a void is identified by 
the subsequent merger of neighboring holes with user-defined criteria (i.e., a percentage of 
volume overwraps with the neighboring holes and a minimum size of hole, ^), where, we use 
^ = I.5R3. As a result, voids identified by the algorithm of FN09 have aspherical shape with 
a finite number of void galaxies within. For each identified void, FN09 found a volume, a 
size and an ellipticity of the best-fit ellipsoid by calculating the inertia tensor from randomly 
scattered particles within the boundary of the void. 

Using a void-finder of FN09, we identify 569, 327, and 804 voids for z = 0, 1, and 2, 
respectively. Here, as we noted earlier, since we are interested in a shape of a tidal field 
near the density minima, we have modified the void finding algorithm of FN09 so that an 
ellipticity of a void can be evaluated from the distributions of galaxies inside a void (i.e., we 
directly calculate an inertia tensor from void galaxies, not from randomly scattered particles 
as in FN09). 

With this modification, the extracted ellipticity PDF has a closer contact with the 
analytic ellipticity PDF derived from the tidal field of the local minima of the density field. 

We measure the linear galaxy biases by comparing the galaxy power spectra of input 
galaxy distributions against linear matter power spectra from CAMB at each redshift. In 
figure 13, we show the linear matter power spectra, Pm{k), from CAMB with the cosmological 
parameters same as the Millennium simulation (i.e., f^m = 0.25, Q\ = 0.75, n^ = 1 and 
(Tg = 0.9) (solid), linear galaxy spectra, Pg{k) = b\Pm{k), (dotted) and power spectra of the 
subset of galaxies from the Millennium simulation (crosses with error bars) at 2; = (top) , 1 
(middle) and 2 (bottom). With fittings up to k^ax = 0.03/i Mpc~ (i.e., using the first two 
data points of Pg{k)), we find bi = 1.0, 1.6, and 2.0 for 2: = 0, 1, and 2, respectively. Here, 
non-linearity of the power spectrum is rather small at this scale, k < 0.03/i Mpc~ . 

In figure 14, we show the ellipticity PDF extracted from the Millennium simulation at 
z = (top), 1 (middle) and 2 (bottom) together with the analytic ellipticity PDF corrected 
for the linear galaxy biases and response functions (red lines). As references, we also show 
the analytic ellipticity PDF (black lines) with (dotted black lines) and without the linear 
galaxy bias (solid black lines). We also show the 1-a range of the histograms calculated with 
bootstrapping method (grey). 

Here, for each redshift, we demand voids have at least 10 void galaxies. Table 2 shows 
redshifts, z, the numbers of voids used to draw histograms, A'v, the minimum numbers of void 
galaxies per void, A'vg,!!!!!!, the median numbers of void galaxies per void, A^vg,medi the linear 
galaxy biases, bz,, the mean density contrasts of void, 6^, the mean Eulerian radii of voids, 
Re = (PiPiPs)^ , the mean ellipticities of extracted voids, e° , and those of the analytic 
PDF, gth°°''y^ and the errors in the mean, Ae"*^^. 

We find from the figure 14 that the analytic ellipticity PDFs convolved with the response 
function (red lines) show a good agreement with the extracted ellipticity PDFs from the 
simulations, with the chi square values of x^ = 23.0, 16.9 and 20.1 with the degrees of 
freedom of 15, 16 and 13 for redshifts of 2; = 0, 1.078 and 2.070, respectively. We test 
the null hypothesis that the extracted ellipticity PDF is drawn from the analytic ellipticity 
PDF convolved with the response function, calculating the probability-to-exceed (PTE) the 
measured x^; 

a = Pn{x)dx, (6-1) 
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Figure 13. We show the hnear matter power spectra, P„i(fc), from CAMB with tlie cosmological 
parameters same as the Millennium simulation (i.e., flm = 0.25, f2A = 0.75, rig = 1 and erg = 0.9) 
(solid), linear galaxy spectra, Pg{k) = b\Pm{k), (dotted) and power spectra of the subset of galaxies 
from the Millennium simulation (crosses with error bars) at z == (top), 1 (middle) and 2 (bottom). 



where Pn{x) is a PDF of x^-distribution with a n-degree of freedom, 



Pn{x) 



X 



n/2-1 -x'^/2 



2"/2r(n/2) 



(6.2) 



We find that our analytic model is consistent with the void ellipticity PDF extracted from 
the Millennium simulation with PTE = 0.12, 0.39 and 0.094 for redshifts of z = 0, 1.078 and 
2.070, respectively. 

Although we have good fits for all redshifts, the extracted ellipticity PDFs are entirely 
dominated by the Poisson noise due to the small number of void galaxies, A^vg ^ 100. To 
see this, we have performed a null test with a top-hat PDF with the mean and the width 
equal to the analytic ellipticity PDF. After convolving the top-hat PDF with the response 
functions, we find an equally good fit with a similar x^- We repeat the null tests for different 
values of the mean and the width of the analytic ellipticity PDF. We find that the extracted 
ellipticity PDF is well fitted with an arbitrary shape of PDF, whose mean ellipticity is the 
same as the one predicted from the analytic ellipticity PDF, convolved with the response 
function. 

6.1 Comments on the Void-Finder Dependence 

In this work, we have assumed that voids are the local minima of density field with three 
dimensionally expanding volume. In the Lagrangian picture, this is equivalent to having all 
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Figure 14. Wc show the eUipticity PDF extracted from the Millennium simulation at z = (top), 
1 (middle) and 2 (bottom) together with the analytic eUipticity PDF corrected for galaxy biases and 
convolved with the response functions (red lines). As references, we also show the analytic eUipticity 
PDF (black lines) with (dotted black lines) and without galaxy bias (solid black lines). We also show 
the 1 — tr range of the histogram, a^^ , estimated via the bootstrap method (gray histograms) . For each 
redshift, the chi-square values are ^ = 21.5, 17.0 and 20.0, and the degree of freedoms are 15, 16 and 
13 for z = 0, 1 and 2, respectively. 
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0.360 


0.008 
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110 
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1.6 


-0.59 


21.9 


0.377 


0.356 


0.012 


2.070 


163 
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-0.60 


19.9 
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Table 2. We show the statistics of voids used to make the histogram in the figure 14. From left to 
right, we have redshifts, z, numbers of voids used to draw histograms, A^v, minimum numbers of void 
galaxies per void, 7Vvg,min, median numbers of void galaxies per void, Avg.mod, linear galaxy biases, 
6l, mean density contrasts of void, 5v, rnean Eulerian radii of voids. Re = ipiViVi)^ (/i~^Mpc), 
mean ellipticities of extracted voids, e°'°^ ^ and those of analytic PDF, g'hoory^ ^^^^ errors in the mean. 



the eigenvalues of the tidal tensor, Tij = 9j(9j^(q), being negative (i.e., A3 > A2 > Ai < 0, 

and Y!i=\\ = ^m)- 

However, as Lavaux & Wandelt [19] pointed out, in a saddle-like density distribution, 
voids can contract along one of the axis, while it expands along the rest of the axes. Lavaux & 
Wandelt [19] studied these spurious voids in the Lagrangian picture together with their own 
void-finding algorithm (DIVA). They found that among those spurious voids, pancake voids 
(i.e., two dimensionally expanding, one dimensionally contracting voids) are as abundant as 
three dimensionally expanding genuine voids, and the shape of the eUipticity distribution 
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function of PLOT needs to be modified. Nevertheless, we focused our work on the three 
dimensionally expanding voids, as the result of PLOT suggests that the void finding algorithm 
of Hoyle & Vogeley [12] mostly finds three dimensionally expanding genuine voids. There 
seems to be an algorithm dependent selection of the type of voids. 

FN09 reported a slight preference for prolate ellipsoids in the ellipticity PDF drawn from 
voids detected in the SDSS DR5, while the analytic ellipticity PDF of PLOT shows a clear 
preference for oblate ellipsoids. FN09 used the best-fit ellipsoid of Jang-Condell & Hernquist 
[15] to estimate the ellipticity of a void with a Monte-Carlo simulation. For each realization, 
they calculated an inertia tensor for each void from test particles randomly spread inside the 
void. In this way, their measurements of ellipticity tend to be more sensitive to the shape 
of the void boundary (i.e., lij oc PiPj). On the other hand, PL07 uses field galaxies inside 
the void to calculate the inertia tensor. An ellipticity calculated from void galaxies is more 
sensitive to the shape of the tidal field around a density minimum, and therefore, has a closer 
contact with the analytic ellipticity PDF of PLOT. Indeed, the mean effective volume of 
voids, Vcs = -f-piP2P3, extracted from the Millennium simulation with the method of PLOT 
is ~ 4 times smaller than that of FN09 (i.e., most void galaxies are within ~ 60% of radius 
from the void center). 

7 Discussion and Conclusion 

We have derived the ellipticity PDF of voids in the redshift space. We have found that 
the redshift space distortion on the shape of voids statistically increases the ellipticities of 
voids, and leaves a prominent feature on the ellipticity PDF as a substantial reduction in the 
probability of having voids with small ellipticity, e < — 4- . This characteristic cutoff in the 
ellipticity PDF can be used as a probe of the growth rate, /(a) = — '^ ina ' o^^e the radial 
density profile of voids is better understood. Acquaviva & Gawiser [1] proposed a model 
independent test of GR by checking the scale dependence/independence of the growth rate, 
f{k,z) from a galaxy power spectrum at small scale, Pg{ks,z), and large scale Pg{kL,z). In 
principle, we can use the redshift space ellipticity PDF of a void to measure the growth rate 
at different scales by binning the observed voids into small and large sizes. 

However, as we have shown in § 5, the biggest limiting factor for the use of ellipticity 
PDF as a probe of cosmology lies in the Poisson noise from a small number of A'vg inside 
a given void. We have found that from a sample of galaxies with a minimum halo mass 
of 2.6 X lO^^/i-^M© for z = and 1, and 8.6 x W^^h'^MQ for 2, we have A^,g ~ 20 void 
galaxies per each identified void. This small number of void galaxies per each void creates 
a significant contamination of the resulting ellipticity PDF so that the shape of the original 
PDF is almost washed-out. 

As we have seen, the biggest limiting factor for the use of ellipticity PDF as a probe 
of cosmology lies in the Poisson noise from a small number of A^vg inside a given void. We 
have found that from a sample of galaxies with a minimum halo mass of 2.6 x 10^'^h~^MQ for 
z = and 1, and 8.6 x lO^^h^^AlQ for 2, we have A^vg ~ 20 void galaxies per each identified 
void. This small number of void galaxies per each void creates a significant contamination 
of the resulting ellipticity PDF so that the shape of the original PDF is almost washed-out. 

Nevertheless, there is a way to overcome the Poisson noise. Recently, Lavaux & Wandelt 
[20] proposed the stacking analysis of the void ellipticity in the redshift space as a way to put 
a constraint on the cosmological parameters via the Alcock Paczynski (AP) test [2]. The AP 
test uses the deviation from the isotropy and homogeneity of the observed structure, which 



28 



is known to be spherical (e.g., galaxy distributions at large scale), to determine D^H. For 
example, when one calibrates comoving distances perpendicular, r_L, and parallel to the line 
of sight, rn, from the angular, 9, and redshift, z, distribution of galaxies with 

. - ^^(^) (7.1) 



{i + z)DA{zy 

c 

one need to assume a reference cosmological model to calculate -D^ (z) and H (z). If the 
observed structure has r± = rn, and we have observed values of the angular, 9, and the 
redshift, z, distributions of galaxies, a combination of D^ and H can be constrained as 
follows, 

(see, e.g., Kazin et al. [17], Shoji et al. [37] for the use of AP test on the two-dimensional 
power spectrum, P{k,fi)). 

In real space, the stacked void with sufficient number of stacked void galaxies should 
have zero ellipticity, while the deviation from zero ellipticity in real space indicates the 
deviation of the assumed cosmology from the true cosmology. In redshift space, since the 
void is elongated toward the line of sight, the stacked void has non-zero ellipticity, which 
can be a tell-tale of the logarithmic growth rate, /(a) = — ^ina Although some useful 
information of void ellipticity will be lost by stacking, in this way, we can see the effect of 
redshift space distortion as a source of anisotropy in the stacked void ellipticity. We think 
that the stacking analysis of the voids in redshift space is potentially a powerful tool to probe 
the cosmology. 
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A Calculation of k(x) 

A mapping between the real-space to the redshift space is 

xl = x\ (A.l) 

^2 ~ ^2 (-^-2) 

xl = xl-fM-±l, (A.3) 

where / = ^i^a ^® ^^^ growth rate, and we set the unit vector, Xg, to be along the line of 
sight direction. Here, u = ^, where v is the peculiar velocity of a void galaxy, and T-L{z) is 
a comoving Hubble rate. 
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We define k(x^) to quantify the strength of the redshift space distortion as fohows 
and we re-write eq. (A. 3) as 



4 = xS(l + «(x^)). (A.5) 

Here, we use a spherical coordinate such that 

xi = xsm6xCos(l)x (A. 6) 

X2 = X sin ^2; sin (/ij; (A. 7) 

X3 = X cos 6x (A. 8) 



and 



ki = A; sin ^fc cos (/>fc (A. 9) 

/c2 = /c sin 0fc sin (/)fc (A. 10) 

/c3 = kcosOk, (A-.ll) 

where x = |x| = y^x^ + Xj + x|, /c = |k| = y^fc^ + k'2 + k'^, and and are the azimuthal 
and inchnation angles respectively. 

/xg • u(xO /xVVcI>(xO 



I A-KGa?pxl 

ik • x^$(k)e* 



/ f d k ,^_ s.r^f^.\J\L■^'^ 



AttGo^pxI J (27r) 
f /""^ r7t /■^'r /"i 



^3 -'0 V^^j JO J-1 



(A.12) 



where p^ = cos^^. Here, k • x = /ex cos 7, where 

cos 7 = COS^fc cos 0a; + COs((/)fc — (px) siu 0fc siu^a;. (A. 13) 

A.l Spherical Void 

As for the spherical void of radius R with homogeneous top-hat density contrast, (5(x), we 
have 

y Jo «;^ 

Therefore, eq. (A.12) becomes 
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We first integrate the angular portion of the integration, J dilfc/i^e*'''''^, expanding the expo- 
nent into Legendre polynomials, 

oo 

^.k.x'- ^ ^(_i)'(2/ + l)ji{-kx)Pi{cosj), (A.16) 

I 

and using the Spherical harmonics addition theorem (a.k.a., Legendre addition theorem), 

4-7r 



m=—l 
I 



m=—l 

Pi{cos 6 k) Pi {cos 9 x) 



I 



{l-m)l^„ 



+ 2 Y. ^-—^^Pr{cosek)Pricosex)cos[m{4>k - </>.)], 
1 1 + ni). 
m=l 



(A.17) 



where PJ^{x) is an associated Legendre polynomial, and it is related to the unassociated 
Legendre polynomial such that 

pr(x) = (-ir(i-x2r/2__p^(^) 

(_-\\m fjl+m 

^ '^-(i-xr/^:^(x^-iy. (A.18) 



2'/! ' ' dx^+' 



Therefore, 



/ dfife^fee^'^-"'' = 5^(-i)'(2/ + l)ji{-kx) I dflkf^kPiicos^) 



/ m=—l 

-0 



^ J_ _ / .^ \3/2 



= -3ihi-kx) l—j Y^\ex,(l>x) = i7riji{kx)cosex, (A.19) 

where we used Y^{d, (j)) = \/ ^ cos9 and ji{—x) = (— l)'j/(x). 
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Finally, we obtain k for a spherical void as 

vrx^ Jo 

A. 2 Spheroidal Void with azimuthal symmetry 

Here, we consider a spheroidal void with an azimuthal symmetry (i.e., an ellipsoid with two 
of the three principal axes have the same length, a = b ^ c). We again choose a spherical 
coordinate, but with the ellipsoidal volume element being projected on to spherical volume 
element such that 

xi = a rsm9cos(j) (A. 21) 

X2 = a rsinOsiiKJ) (A. 22) 

X3 = c rcos9, (A. 23) 

where < r < 1 and d'^x = a^c r'^dr d(p dcosO. 

First, let us derive the Fourier transform of the spheroidal void with a homogeneous 
top-hat density contrast, (5(x) as follows, 

/■I r 

(A.24) 





S{k) = 6 1 d^xe-'^'' = o?cl f r^dr f dne-'^'' 


where 






k • X = /cr[a sin 0fc sin 0^. cos ((/)fc — (j)x) + c cos 0^ cos 6 


Here, we define 






cos 7a = sin 9k sin 0^ cos((/>fc - (px) 




cos 7c = cos Ok cos 6x , 



(A.25) 



(A.26) 
(A.27) 



and 



cos 7 = cos 7a + cos 7c (A. 28) 

COST] = — cos 7a + cos 7c, (A. 29) 

where cosry is obtained by changing a sign of Ox- In terms of cos 7 and cos 77, we have 

k • X = -kr[{c + a) COS7 + {c — a) cos?/]. (A. 30) 
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Again, we use the spherical harmonics addition theorem (eq. (A. 17)). Note that due to the 
definition of eq. (A. 29) and from eq. (A. 17), we see that Pi{cosr]) = Pi{cosj), and therefore, 



-jkx 



exp[—-ikr{c + a) C0S7] exp[—-ikr{c — a) cost/] 
^(-i)'(2/ + l)ji (hric + a)] Piicosj) 

J^i-if {21' + l)ji, (hr{c-a)) P^cost) 
i'=o ^^ ^ 

1=0 1'=0 ^ 

X Ji' ( -^'^(c - «) ) Pi{cos-f)Pi>{cosj). 



From the Legendre addition theorem, we have 
dQ^Pi (cos 'y)Pi> (cos 7) 

I I' 
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(2/ + l)(2/ 
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'll"J—mm' ; 
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then we have 



d!^^ e-**^-^ 
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i=0 /'=0 
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X Ji' [^f^'i^ic- a)j / dO^.P/(cos7)P//(cos7) 

= IStt^E E (-1^"^^' (^^r(c + a))i, (^^r(c - a)) 
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X Yr{6k,<t^k)Yf^{ekAk). 



(A.31) 



(A.32) 



(A.33) 



33 



Therefore, 



6{k) = a^c6 / r'dr / dOe 

oo I 



-ikx 
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1=0 m=-l 



X / r 



drji -kr{c + o) ji -kr{c - a) 



1 



and integrating over r G [0, 1], we have 



^(k) 



Svr^aJ 



oo I 



\J (? — a'^k'^ 



^ ^ (-i)'+'"yr(^;L.<^.)>^r'"(^fc''^/^) 



/=0 m=~l 



(c + a)J,_i f-A;(c + a)j J^^i l-k{c-a) 
- (c-a)J;_i f-/c(c-a)j J;_^i (-k{c + a) 
We Taylor expand the above equation for 6(k) about e ~ up to a linear order. 



5(1:) -.-^^"V (^-y^'^' 

~ ^' ^ 22'-3/2r(/ + 3/2) 



kaeJi^i{ka) + {2 + (/ - l)e} J^_^3 (A;a) 



For a prolate ellipsoid, where a < c and e = 1 — |, 

xn ^ /I X 3 [ ii(M , f jijka) . ., ^ 



j2(A:a) 



fca \ ka 

For an oblate ellipsoid, where a > c and e = 1 — |, 

' ji{ka) ( ji{ka) 



(5(k) ~ 47r(^a'^ 



ka 



ka 



+ 0{e^ 



+ 0{e'). 



Therefore, eq. (A. 12) becomes 



k;(x) 



A-Ka^f5 



dk 



X3 Jo (2vr)' 
dUkCosOke''^''^, 



rik 



ji{ka) f jijka) . 



ka 



ka 



where 



dQu cos 



(A.34) 



(A.35) 



(A.36) 



(A.37) 



(A.38) 



(A.39) 



OO OO r~. / 1 \ 

0^,.k.x = ^ ^ ^^(_,)W'(2/ + l)(2/' + l)ji (--kr{c + a) 
1=0 i'=o ^ ^ 

X ji' i-2^r{c-a)] / (iJ]fcy]°(6'fc, (^:>fc)-PKcos7)-P^'(cos7)• 



(A.40) 
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Using the addition theorem given in eq. (A. 17), we rewrite the above equation, 



dilfc cos Oke"^'^ 



oo oo i I' r: y ^ 



/=0 l'=Om=~lin'=~l' 



1 



\m+m' \/ — m./ 



X ji' l--kr{c-a)\ (-ir+"^yr"(^.,<Ax)V"(^x''/'-) 

X / dnkY^^iOk, (Pk)Yr{Ok, 4>k)Yp'{ek, 4>k). 



(A.41) 



Here, 
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-kr{c + a) j iv ( --kr{c - a) 



-21' 



27r 



(-1) 



l+l' 



r(^' + i: 



-{kra) 



I'-- I 



{2 + (/ + 2l')e}Jij^i{kra) - kraeJi_^3{kra) 
l+V 



^-2(r+i) V J ^j^^^Y^V [|2 + (/ + 2l')e}ji{kra) - kraeJi+,{kTa)\ 



(A.42) 



for a prolate elHpsoid {c> a), and 



-kr{c + a) j jV f --fcr(c - a) 



r(/' + |)^ ' 

2a'+i) \/vr(-l)^ 



(2 — le)Jij^i{kra) + kraeJi_^3{kra) 



2-^^' +.;^_v^(^^^y'/ [(2 - k)jiikra) + A:raeji+i(A;ra)] 



T{1' + 



(A.43) 



for an oblate ellipsoid {c < a). In order to truncate the equation at linear order in e, we have 
I' < 1, such that 



Ji I -^kr{c + a)j ji, (--kr{c-a) 

i(-l)' [(2 + le)ji{kra) - kraeji+i{kra)] 
-\{-l)^kraeji{kra) 

for a prolate ellipsoid, and 

31 i-^kr{c + a)j ji> {--kr{c - a] 

i(-l)' [(2 - le)ji{kra) + kraeji+i{kra)] 
~ ■^ ^{ — iykraeji{kra) 



for /' = 

for /' = 1 

otherwise 



for I' = 
for I' = 1 
otherwise 



(A.44) 



(A.45) 
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for an oblate ellipsoid. 

As for the integration over the solid angle, ftk, we have 



dnkY,\9k, MYnOk, MYP {Ok, (t>k) 



3(2/ + l)(2/^ + l) / 1 / /' \ / 1 / /' 

Att loOO/lomm''' 



(A.46) 



where the Wigner-3J symbol must satisfy the triangular inequality such that |1— /| < /' < 1+^, 
and + m = —m! in order to have a non-zero value. 

As for eq. (A. 41), since I' < 1, we have a limited number of non-zero terms. We have 
{I, V, m, m') = (0, 1, 0, 0), (1, 0, 0, 0), (2, 1, -1, 1), (2, 1, 0, 0) and (2, 1, 1, -1) as the non-zero 
terms. We rewrite eq. (A. 41) with a finite number of terms truncating 0{e^) and higher as 
follows. 



/ 



dQ.kCosOke''^'^ ~ 



^/3 



±hraeJo{kra)YS{e^, (t>^)Y^{e.Ax) 



+ - {(2 ± e)ji{kra) ^ kraeJ2{kra)} Y^{e,,^,)Yl>{e,,cb,) 



zb —kraeJ2{kra) 



■^y; 



'xAx)Y2\e^,(l)^) 



+ 



\j\Yl{e^,(t>^)Y^\e^Ax) - ^Y^{e^,^.,)Y^{e^,ci)., 



(A.47) 



Here, in the limit of e — )■ 0, we recover 

\2 



d^i,cos9ke''^''^ 



^-^iji{kra)Y^{e^, (t):,)Y^{e^, (j)^) = Anijiikra) cos ( 



(A.48) 



as in eq. (A. 19). 
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Finally, from eq. (A. 39) and (A. 47), and truncating ©(e^) and higher terms, we have 



2^3 Jo V'^Y 



ik 



jijka) 
ka 



+ e 



Jijka) 
ka 



J2{ka) 



dUk cos ^fce' 
8f6a^ 



jkx 



VSx-i 



-eY^{9,, ^^)Y^{9.„ </>,.) / dk kji{ka)jo{kra) 
6 Jo 



+ 



2 + e 
2a 



Yo'ie.. 



X, •i'x)Yl {dx,4'x) 



dk ji{ka)ji{kra) 



-eY,\e,,^,)Y,°{ 



,) / dk kji{ka)J2{kra) 
Jo 



+ 



h 1 Vf ^i~'^^"' '^-)^2'(^x, <^x) + y|n'( 



'X 1 (px ) ^2 



- -Tf'ii^x, <Px)Y^{Ox, <Px) I / dk kjiika)J2{kra) 
+ -Y^{0^, (^x)Y^^{ex, 4>x) / dk h{ka)3i{kra) 



eY^{e^,ck,)Y^{e^,, 



dk kJ2{ka)ji{kra) 



Taking integrations over k (see appendix §B), 

f6 ar cos 9^ 



K X 



6 X3 

-f(l-2.) 



(1-e) 



Similarly, for an oblate spheroid, 



fl 



(A.49) 



(A.50) 



(A.51) 



As we compare against numerical calculations of k(x), we see k(x) ~ — 4-(l =p e) fits the 



III 



K(x)~-^(l + 2e) 

ms 
result better for both prolate and oblate spheroids. 

We see how real space geometry of a given void affects the peculiar velocity field within 
the void by comparing Eqs.(A.20) and (A.50). As for a prolate spherical void, we set the 
longest principal axis, a < c= a/(l — e), to lie along the line of sight. As discussed by Icke 
[14], originally aspherical voids tend to become a sphere with a slower expansion speed along 
the longest principal axis, and vice versa. For a given prolate spheroidal void with its longest 
principal axis along the line of sight, peculiar velocity along the largest axis become smaller 
as we increase ellipticity. Since the deformation of the void shape in the redshift space, k(x), 
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is proportional to the peculiar velocity, an aspherical void in real space with the longest axis 
being along the line of sight experiences less deformation compared to an originally spherical 
void. 

Also, we see that the degree of deformation of an aspherical void in the redshift space, 
k(x), is independent of the size of the void, but dependent only on the linear growth factor, 
/ = ^^° ^ , the mean density contrast of the void, 6, and the ellipticity of the void in the real 
space. 

B List of y^™, ji and Wigner-3J symbols used in appendix A 

Here, we collect the properties of the spherical harmonics, y^'"(0, 0), spherical Bessel func- 

I , that are 
used in our calculations. 



20F 



z V -^vr 



Ji{z) 



smz 

z 
sin z cos z 



z^ z 



B.l 
B.2 
B.3 
B.4 

B.5 
B.6 
B.7 

B.8 
B.9 

r(3/2) = ^ (B.io: 

r(5/2) = ^ (B.ll 

r(7/2) = 1^ (B.12 

r(9/2) = 1^ (B.13 



1\ 3 

32KZ) = [ ^--jsmz-^cosz 



r(i/2) = V^ 
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1 1\ i_ 

OOOJ " V3 

1 1 0\ 1_ 

OOOy' " V3 

1 2 1\ _ /T 

oy V 15 

1 2 1\ 1 



o-iiy Vw 

1 2 1 \ 1 



For < r < 1 and a > 0, 



1 -ly ^ 



vrr 
dk ji{ka)ji{kra) = -— 
ba 





°° vrr 

dfc k J2{ka)ji{kra) = —^ 

j dk k ji{ka)jo{kra) = -^ 

dk k ji{ka)J2{kra) = 
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B.15) 
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B.18) 

B.19) 
B.20) 
B.21) 
B.22) 
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